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Abstract. In the perturbative QCD with N c — > oo equations for the amplitude of the 
nucleus-nucleus scattering are derived by the effective field method. The asymptotic form of 
the solution is discussed. It is argued that in the high-energy limit the total nucleus-nucleus 

m : 

cross-sections become constant and purely geometrical. 
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1 Introduction 

^p Hi In the framework of the colour dipole model of A.H.Mueller [1,2] it follows that in the high- 

colour limit N c — ► oo the scattering on a heavy nucleus is exactly described by the sum of fan 
diagrams constructed of BFKL pomerons, each of them splitting into two [3]. The equation 
for the sum of BFKL fan diagrams was first written by I.Balitsky [4] in his original operator 
expansion formalism. Then it was rederived by Yu.Kovchegov [5] in the colour dipole frame- 
work and by the author by directly summing the BFKL fan diagrams [6]. The perturbative 
solution of this equation in the region of small non-linearity (outside the saturation region) 
was studied in [7]. Asymptotic estimates of the solution were presented in [8]. Finally in 
[6] the exact solution of the equation was obtained by direct numerical methods. The main 
physical results following from these studies are that, first, at high rapidities Y the hA total 
cross-section saturates to its geometrical limit 2ttR\ and, second, the gluon density in the 
nucleus aquires a form of a soliton in Y — Ink space moving towards higher momenta with 
nearly a constant velocity as Y increases. This last property supports the applicability of the 
perturbative treatment, since the well-know diffusion of the BFFKL pomeron towards small 
momenta results to be stopped. 

In the present paper we attempt to generalize these results to nucleus-nucleus (AB) scat- 
tering. In this case, in the N c — > oo limit the total amplitude is given by the sum of all tree 
diagrams constructed of BFKL pomerons and the triple pomeron vertex. In contrast to the 
hA case, the vertex now describes not only splitting of a pomeron into two but also fusion 
of two pomerons into one. The diagrams for the amplitude accordingly result much more 
complicated than the fan diagrams relevant for the hA case. However, using the effective field 
theory methods developed for summing such diagrams long ago [9,10], one can construct an 
equation which describes the AB amplitude. Naturally this equation (in fact a pair of equa- 
tions) results much more complicated than for the case of hA scattering. Its exact (numerical) 
solution does not seem realistic. However simple asymptotic estimates, analogous to the ones 
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made in [7,8 ], show that the total AB cross-section tends to its geometrical limit at high 
rapidities similar to the hA case. 

Unfortunately the gluon density in the overlapping area cannot be found from these 
estimates, but rather requires knowledge of the solution in more detail. We leave this problem 
for future studies. 



2 Effective field theory for AB scattering 

At fixed overall impact parameter b the AB amplitude A(Y, b) can be presented as an expo- 
nential of its connected part: 

A(Y,b) = 2is(l - e ~ T(Y ' b) ) (1) 

The dimensionless T is an integral over two impact parameters b A and b B of the collision 
point relative to the centers of the nuclei A and B: 

T(Y, b) = J d 2 b A d 2 b B 5 2 {b -b A + b B )T(Y, b A , b B ) (2) 

As mentioned in the Introduction, in the perturbative QCD with N c — > oo the amplitude 
—T(Y, b A , b B ) is given by a sum of all connected tree diagrams constructed of BFKL pomerons 
and the triple pomeron vertex. More concretely, in these diagrams a line ("propagator") 
connecting two points j/i,ri and y 2 ,r 2 corresponds to one half of the forward BFKL Green 
function [11]: 

G(y 1 -y 2 ,r 1 ,r 2 ) = — 2 e(y 1 -y 2 )j:j *» J < < (3) 

where 4>i t2 and are the azimuthal angles and 

w (!,) = ^(^(l)-Re^(l/2 + M/)) (4) 

7T 

are the BFKL levels. Due to the azimuthal symmetry of the projectile and target colour 
densities one may retain only the term with zero orbital momenta n = in (3). 

The interaction between the pomerons is realized via the triple pomeron vertex. It is non- 
local and not symmetric in the incoming and outgoing pomerons. Its form for the splitting 
of a pomeron into two was established in [2,12,13]. At N c — > 00 for the transition 1^2 + 3 
the three BFKL Green functions are connected by it as follows (see Fig. la) 



vr J r^r^r^ 



Here it is assumed that the operator Vi acts on the left. The form of the vertex for the fusion 
of two pomerons into one is actually not known. However, the symmetry between target and 
projectile prompts us to assume that for the inverse process 2 + 3^1 the BFKL functions 
are to be joined as (Fig. lb) 

4a 2 s N c f d 2 rid 2 r 2 d 2 r^ 2 , , , , 4 4 , 
/ 2-2-2 6 (. r i+ r 2+r3)G(y 2 -y,r 2 ,r 2 )G{y 3 -y,r 3 ,r 3 )V 1 r 1 G{y-y 1 ,r 1 ,r 1 ) (6) 



Finally we have to describe the interaction of the pomerons with the two nuclei. The 
BFKL Green functions corresponding to the external legs of the diagrams are to be integrated 
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with the colour density of each nucleus. We take the target nucleus at rest, that is, at rapidity- 
zero. Then each outgoing external BFKL Green function is to be transformed into 

g 2 AT A (b A ) J d 2 r'G(y, r, r')p N (r') = J dy'd 2 r'G{y - y', r, r')r A (y f , r') (7) 

where pat is the colour density of the nucleon, T A is the profile function of the nucleus A and 
we define 

r A (y, r) = g 2 AT A {b A )p N (r)5{y) (8) 

(with dependence on b A implicit). Similarly each ingoing BFKL external Green function is 
transformed into 

J dy'd 2 r'G(y' - y, r', r)T B (y', r') (9) 

where 

tb(v, r) = g 2 BT B (b B ) PN (r)5(y - Y) (10) 

To find the amplitude one has to sum over all connected diagrams with M ingoing and N 
outgoing lines, corresponding to M interactions with the projectile and N interactions with 
the target, divided by MINI. 

It is trivial to see that this sum exactly corresponds to the sum of tree diagrams generated 
by an effective quantum theory of two pomeronic fields <E>(y, r) and 3>'(y, r) with action 

S = S + S I + S E (11) 

consisting of three terms, which correspond to free pomerons, their mutual interaction and 
their interaction with external sourses (nuclei) respectively. 
To give the correct propagators Sq has to be chosen as 

S = 2 J dy 1 d 2 r 1 dy 2 d 2 r 2 <5>(y 1 ,r 1 )G~ 1 (y 1 - y 2 ,r 1 ,r 2 )&(y 2 ,r 2 ) = 2($\G- l \&) (12) 

where (|) means the integration over y,r. Note that the sign of Sq corresponds to the 
substitution of the conventionally defined field variables $ and 

$t_>j$t (13) 

which allows to make all terms of the action real. 

According to (5), (6) the interaction term Si is local in rapidity 

Si = ^Nc f dy I ^!l|!r^ ^ (ri+r2+r3) r v 4 r 4 $ t (y;ri)$( y ,r 2 )$(y,r 3 ) + c.c] (14) 

The overall sign combines the initial factor i and i 3 from the substitution (13). 
Finally the interaction with the nuclei is local both in rapidity and coordinates: 

Se = ~ J dyd 2 r\^(y,r)T A (y 1 r) + $ f (y, r)r B (y, r)] (15) 

The minus sign comes from the initial i and the substitution (13). 

The amplitude T(Y, b Al b B ) is then expressed through a functional integral 

Z = j D<5>D&e s/ ^ (16) 
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where jjl is an arbitrary mass scale necessary to adjust the dimensions of various parts of the 
action. Keeping only the connected diagrams one finds 

T(Y,b A ,b b ) = -^ 2 ln^- (17) 

where Zq is the value of Z for Se = 0. Functional integral Z is to be calculated in the 
classical approximation to retain only the tree diagrams. This gives 

T(Y,b A ,b B ) = -S E {*,&} = J d 2 r[<S>(0,r)f A (r) + &(Y,r)t B (rj\ (18) 

where $ and $^ are the solutions of the classical equation of motion and f's are (8) and (10) 
with the S functions of rapidity dropped. We see that the result is independent of the scale 
/j,, as it should be. 

3 Equations of motion 

Before writing out the classical equation of motion, we transform the functional integral (16) 
to new variables in which the non-locality of the Lagrangian becomes substantially reduced. 
We put 

$(y, r) = r 2 (j)(y, r), $ f (y, r) = rV(r, y) (19) 
In these variables the interaction term becomes 

S i = ^^ £ / d y J d 2 r 1 d 2 r 2 d 2 r :i 5 2 (r 1 +r 2 + r 3 )[K 1 ^(y,r 1 )(t } (y,r 2 )cl)(y,r 3 ) + c.c] (20) 

where (dimensionles) operator K has the form 

K = r 2 Vjr 2 (21) 
Now we transform (20) to the momentum space to obtain 

Si = 40^ j dy j ^^ ]{y ^ q) ^ q) ^ q) + cc ] (22) 

with operator K local in the momentum space 

K = V 2 qq 4 V 2 q (23) 

As we observe, the interaction has become local in the momentum space. 
Now we turn to the free part So- In new variables it takes the form 

S = 2(<j ) \r 2 G- 1 r 2 \^) (24) 

It was shown in [6] that 

r 2 y\G(y,r u r 2 ) = g{y,ri,r 2 )rl (25) 

where g is the Green function of the BFKL equation for the so-called semi-amputated wave 
function: 

' } +H 1 )g(y,r 1 ,r 2 )=5(y)6 2 (r 1 -r 2 ) (26) 



Here Hi is the BFKL Hamiltonian [11] acting on r±. We rewrite (26) in the operatorial form 

Kr- 2 Gr- 2 =(Jj + H) ' (27) 
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wherefrom 

^G'^K-^-^ + H (28) 
dy 

and finally 

r 2 G~ 1 r 2 = (^ + H\K (29) 



dy 

Since both G(y, r±,r2) and g(y, r\,r-2) are symmetric in ri,V2 we also find in the same manner 

r 2 G- l r 2 = k(^- + h] (30) 



dy 

so that K commutes with H. Using (29) and (30) in (24) we see that the free part has become 
local in rapidity and can be expressed via the BFKL Hamiltonian H: 

S = 2(<I>\K(J- + Hj\^) (31) 

This part remains non-local both in the coordinate and momentum spaces due to the non- 
locality of H. 

The interaction with the nucleus part in the new variables takes the form 

S E = -(w A \0)-(^\w B ) (32) 

where, in the coordinate space, 

wa,b(v, r) = r 2 T A)B (y, r) (33) 

Now, with the action nearly completely local (except for the Hamiltonian term) we can 
write out the equation of motion. We find 

5S 2K (J- + Hj tf(y, q ) + (K^ 2 (y, q) + 2</>(y, q)K^(y.q)) - w A (y, q) = 

(34) 



and 
5S 



~^ = 2K (-|- + Hj <f>(y, q) + (K0 2 (y, q) + 2^(y, q)K<f>(y.qj) - w B (y, q) = 

(35) 



S<f>Hy 

Applying operator (1/2)K~ 1 from the left we find our final equations of motion 

d „\ , + , N 2a 2 N c / , +2 



— + Hj <j>\y, q) + (^(y, q) + 2R- 1 [cP(y, q)K<j>\y.q)]) - l -K- l w A {y, q) = (36) 



and 

~ + Hj </>(y, q) + (0 2 (y, q) + 2K~ 1 [<f>\y, q)K</>(y.q)]) - ^K^wsiy, ?) = (37) 

As we see the resulting equations are rather complicated, since they involve nonlocal 
terms, bilinear in <p and $ , which interconnect the two equations. Summing fan diagrams 
in the hA case in fact leads to the the same equations, in which however wb = 0. Then one 
immediately finds that = identically, which converts the first equation into 

(j^ + Hj tf{y, q) + ^V 2 (y) q) ~ ^K-'wAiy, q) = (38) 
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This local equation is just the one studied in [4-6]. 

To conclude this section we present the non-local terms in Eqs. (36) and (37) in a more 
explicit way. To this end we calculate the kernel of the operator K^ 1 in the momentum 
space. We have in the coordinate space 

R- 1 = r- 2 V; 4 r- 2 (39) 

Using the identity [6] 

ViG(0,r 1 ,r 2 ) = 8 2 (r 1 -r 2 ) (40) 
we can write the kernel of K^ 1 in the coordinate space as 

K- l { ri ,r 2 ) = r^ 2 G(0,n,r 2 )r 2 2 (41) 

Fourier transforming (41) and using (3) we find the kernel in the momentum space: 

K-\qi,Q2) = J ^^e^-^G(0,n,r 2 ) = \ J ^-j^ Jfo qi )I* („, q 2 ) (42) 

where 

I(u, q) = drr-^Mqr) = 2-"V 1 +*» ^ " iv) (43) 
JO 1 (1/2 + iv) 

Doing the integral over v we finally find 

where ?>(<) = max(min){gi, ri- 
nsing (44) we can rewrite the nonlocal term in Eq. (36) as 

±o? s N c j (ln^ + l) Hy,QiMiM^(y,Qi) (45) 

where (?>(<) = max(min){(7, qi}. The nonlocal term in Eq. (37) is obtained by complex 
conjugation. 

4 The total cross-section 

The obtained equations which determine the classical fields <f> and <ft are very difficult to 
solve even numerically. The trouble lies not so in the their explicit non-locality, but in the 
appearence of two different sources at two different rapidities. Due to conditions (f) = for 
y > Y and <f>^ = for y < and the (5-like dependence of the sources on rapidities, one can 
drop the sources in Eq. (36) and (37) substituting them by conditions 

4>\y, q) y =o = K~ l w A {q), 4>{y, q) y =Y = K~ l w B {q) (46) 

In the hA case one has only the first of these conditions, which converts Eq. (36) into an 
evolution equation in rapidity, relatively easily solved by conventional methods. As mentioned 
the non-local term is zero in this case, but its presence would only slightly complicate the 
solution. After all the BFKL Hamiltonian is also non-local (although linear). 

For the nucleus-nucleus scattering we have to solve homogeneous Eqs (36) and (37) with 
both conditions (46) imposed upon the solution. The Cauchy problem is thus transformed 
into an essentially more difficult boundary problem. A possible method of the solution is to 
transform Eqs. (36) and (37) into a system of two non-linear integral equations in combined 
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rapidity-momentum space, which one may try to solve by iterations. Our experience in the 
hA problem shows that for the solution to have a reasonable precision one requires at least 
800 points in the momentum and 400 points for 5 units of rapidity. Thus to find the amplitude 
for say Y = 15 one has to perform 1200 x 800 3 ~ 10 10 operations per iteration. On top of 
that the convergence properties of the iteration procedure is unknown. 

Here we shall not attempt at solving Eqs. (36) and (37) with any reasonable degree of 
precision at all values of rapidity and momentum. Instead we shall again use our experience 
with the case of hA scattering (fan diagrams), where at any fixed momentum and y — > oo 
the solution <f>*(y, q) aquires a simple form, independent of the target properties 

4>Hy,<l)y^oo = Jln^p (47) 

with lnQ(y) ~ 2.3A(oi s N c /ir)y. Note that (47) is not the solution at all y and q. In particular 
(47) is not valid at q ~ Q, which is just the region which determines the gluon density. How- 
ever (47) is sufficient to establish that the hA scattering cross-section tends to its geometric 
limit at high y [7]. 

Our guess is that also in the nucleus- nucleus case function <f>'(y,q) aquires the form (47) 
at large rapidities y ~ Y and cj>(y, q) aquires the same form with y — > Y — y. To support this 
behaviour we are going to demonstrate that in these limits the non-local terms in Eqs. (36) 
and (37) can be neglected, so that the equations decouple and become similar to the hA case. 

Actually the demonstration is quite simple. Take Eq. (36) at large y and put the con- 
jectured asymptotics (47) into it. According to (45) the mixing non-local term will then be 
given by 

a ' N * I T^T fa - + l) HV, *i) VfrfV? In (48) 
■ (2tt) 2 ^ V q< J q\ 

Note that (f>(y,qi) enters at large y and small Y — y. Its exact form is unknown but we 
can safely assume that it rapidly falls with q\ similar to the inhomogeneous term K~ 1 wb 
from which it is separated by a relatively small distance in rapidity. Action of the operator 
K = V 2 q 4 V 2 on the asymptotic form of <ft however gives zero. In fact 

V 2i n ^l = _ 2 ^ ((?l) 
Qi 

Subsequent integration over q\ gives zero at any finite q due to factor qf. Thus the non-local 
term is zero in Eq. (36) at y ~ Y — > oo. Therefore the equation aquires the same form as for 
hA scattering (which implies that fan diagrams going from top to bottom dominate). This 
means that the asymptotical behaviour (47) is indeed true. The same result is found for the 
non-local term in Eq. (37) at small y and Y — y — ► oo assuming the asymptotic form (47) 
for (f)(Y — y,q). Its meaning is that fan diagrams going from bottom to top dominate in this 
limit. 

Functions $(y, r) and &'(y,r) which actually determine the amplitude according to (18) 
are related to 4>(y,q) and 4>^(y,q) by 

*(v,r) = -/ J^e iqr V 2 q <f>(y,Q) (49) 

and similarly for the conjugated function. Using the asymptotical expression (47) we then 
get at high Y 

$(0,r) = \9(R B -b B ), &(Y,r) = \o(R A -b A ) (50) 
9 9 
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The two 6 functions appear because according to Eqs. (36) and (37) = = 0) when 
wb = (wb = 0), that is for bs > Rb Q>A > Ra)- 

Putting (50) in (18) we obtain the connected part of the amplitude at Y » 1 as 

T(Y, b A , b B ) = AT A (b A )9(R B - b B ) + T B (b B )9(R A - b A ) (51) 
It results independent of Y. After the integration over bs and 6g we find 

T(Y, b) = J d 2 b A d 2 b B 5 2 (b -b A + b B )9{R A - b A )6{R B - b B )[AT A {b A ) + BT B (b B )] (52) 

According to (1) the total AB cross-section is given by 

o to \Y) = 2 J d 2 b(l-e- T W) (53) 

/.From (52) and (53) one concludes that at large Y the cross-section does not depend on Y. It 
saturates at a value which is purely geometrical and for A >> 1 or /and B >> 1. approaches 
the black disk limit in the overlap area. 



5 Conclusions 

We have derived a pair of equations which describe the nucleus-nucleus scattering in the 
perturbative QCD with a large number of colours (or, alternatively in the quasi-classical 
limit, or, in the limit A, B — > oo). The equations contain mixing terms which are both non- 
linear and non-local. In contrast to the hA case the equations are to be solved with given 
boundary conditions at rapidities both of the projectile and target, which complicates the 
solution enormously. 

However the asymptotical form of the solution at fixed momentum and large rapidities is 
shown to be the same as for the hA case. This allows to demonstrate that at large rapidities 
the total AB cross-section becomes independent of energy and given by purely geometric 
considerations. At large A or/and B it corresponds to the scattering on the black disc in the 
overlap area. 

Going to particle production in AB collisions, the situation at central rapidities seems to 
be rather simple. The inclusive cross-section will be described by diagrams like shown in Fig. 
2, with the target and projectile parts joined by a single pomeron, from which the observed 
particle is emitted. Evidently this contribution is just a convolution of the production vertex 
with two gluon densities of the projectile and target, each one corresponding to fan diagrams 
and found in [6]. At rapiditiy distances from the target or projectile 5y ~ 1/A where A is 
the pomeron intercept the problem does not look so simple, since the AGK rules are rather 
complicated in this region (see e. g. [14]) and the usual cancellation of all diagrams except of 
the structure shown in Fig. 2 is not at all obvious. We leave this problem for future studies. 
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7 Figure captions 

Fig. 1. The triple pomeron vertex for the splitting of a pomeron into two (a) and fusion of 
two pomerons into one (b). 

Fig. 2. The generic diagram for the inclusive particle production in the central region in AB 
collisions 
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Fig. 2 



